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Abstract 

Collinear factorized perturbative QCD model predictions are compared for p + Pb at A.AA TeV 
to test nuclear shadowing of parton distribution at the Large Hadron Collider (LHC). The nuclear 
modification factor (NMF), RpPh{v = 0,pT < 20 GeV/c) = dnpp\^/ {Nco\\{b)dnpp), is computed 
with electron-nucleus (e + A) global fit with different nuclear shadow distributions and compared 
to fixed shadow ansatz used in Monte Carlo Heavy Ion Jet Interacting Generator (HIJING) 
type models. Due to rapid DGLAP reduction of shadowing with increasing used in e + ^ global 
fit, our results confirm that no significant initial state suppression is expected {RpPhiPx) = 1±0.1) 
in the pT range 5 to 20 GeV/c. In contrast, the fixed shadowing models assumed in HIJING 
type models predict in the above pr range a sizable suppression, RpPhipr) = 0.6 — 0.7 at mid- 
pseudorapidity that is similar to the color glass condensate (CGC) model predictions. For central 
(-^coii = 12) p + Pb collisions and at forward pseudorapidity (r/ = 6) the HIJING type models 
predict smaller values of nuclear modification factors {RpPhiPr)) than in minimum bias events at 
mid-pseudorapidity {rj = 0). Observation of RpPh{PT = 5 — 20 GeV/c) ~ 0.6 for minimum bias 
p + A collisions would pose a serious difficulty for separating initial from final state interactions in 
Pb+Pb collisions at LHC energies. 

PACS numbers: 12.38.Mli, 24.85. +p, 25.75.-q, 24.10Lx 
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I. INTRODUCTION 



In this paper we compare predictions for moderate pt < 20 GeV/ c observables in p + Pb 
reactions at 4.474 TeV at the Large Hadron Colhder (LHC), that will help to discriminate 
between models of initial conditions assumed in Pb+Pb collisions at 2.76 A TeV. The possi- 
bility that the first data on p + Pb may be taken soon, with a potential high physics payoff 
motivates this paper. All models details are extensively discussed in the literature and we 
focus only on the updated nuclear modification factor, i?pPb(?7,PT, b) = dupPh/ {Nco\\{b)dnpp) , 
predictions testable with a short AAA TeV run. In minimum bias (MB) and central (0—20%) 
p + Pb collisions the average number of binary nucleon-nucleon (NN) interactions (with an 
inelastic cross section crj^^ ~ 65 mb) is ^ 7 and N^^l^ ^12 respectively. 

This control experiment has long been anticipated to play a decisive role in helping to 
deconvolute initial and final state interaction effects in Pb+Pb reactions at the LHC The 
d + Au control experiment at 0.2^4 TeV played a similar critical role for Au+Au at the 
Relativistic Heavy Ion Collider (RHIC) in 2004 jl, 2|. The importance of p + Pb was also 



emphasized in the 2007 Last Call for LHC compilation of predictions [3[ and many other 
works j4-10]. 

The open problem after the first very successful LHC heavy ion run in 2010 ll| remains 
how to deconvolute nuclear modification effects due to initial state and final state effects. 
Without a clear calibration of the magnitude of initial state suppression of the incident 
nuclear partonic flux it is not possible to draw firm conclusions about the properties of the 
quark gluon plasma (QGP) phase of matter produced at the LHC. At RHIC the same prob- 
lem was resolved at mid-rapidity by the observation of no appreciable nuclear modification in 
d + An control experiments in 2003 (see also Figure [3] below) . An approximately factor 
of four suppression of moderate pt mid-rapidity pions observed in Au+Au at RHIC could 
then be interpreted as due to final state jet energy loss in a high opacity QGP produced in 
central Au+Au collisions at RHIC. 

At the LHC the initial flux is much more uncertain than at RHIC because of the higher 
density of partons at an order of magnitude smaller fractional momenta x = 2pt/\^ < 10~^. 
At high initial densities all models predict a breakdown of additivity of the nuclear parton 
distribution functions (nPDF). The magnitude of the breakdown however varies greatly in 
the literature in both coUinear factorized approaches and kx factorized parton saturation 



2 



model approaches 
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10|. Therefore, even a rough first experimental constraint from 



p+Pb interactions would have high impact on the development of nuclear collision modeling. 



II. NUCLEAR SHADOWING AND JET QUENCHING AT LHC ENERGIES 

Nuclear shadowing of quark and gluon nPDFs at large x > 0.01 and moderate is 
well constrained from e + A and lower energy p + A data. Global fit parametrizations of 
the nuclear parton distribution functions (nPDF) are available flSQ. DGLAP evolution 
[l^ to higher predicts a rapid reduction of shadowing effects and therefore only modest 
modifications of RpPhipr) = 1 ± 0.1 for pr > 5 GeV/c have been predicted 

HQ- As 

emphasized in [lO], any observed significant modification from unity would be inconsistent 
with most current nPDFs and therefore pose a severe challenge to conventional coUinear 
factorized QCD approximation to high-po" processes not only in p + Pb but even more so in 
Pb+Pb collisions. We continue here to investigate this central thesis. 

At RHIC there is clear evidence at high rapidities, where small fractional parton momenta 
X ~ 10^'^ similar to central p + Pb are probed, that binary collision scaling of collinear 
factorization breaks down. Color Glass Condensate (CGC) kx factorization models have 



been developed to explain these deviations Q, [Q and nearly identical RpA ~ 0.7 ± 0.1 
nuclear modifications factors were predicted in Ref. [l^ (KKT04) for forward rapidities at 
RHIC and mid-rapidity at the LHC. 

However, collinear factorized approaches with DGLAP evolved nPDF appear to provide 
an alternate explanation of forward single inclusive yields at RHIC fl. At LHC energies 
we can differentiate between these explanations because the collinear factorized approach 
predicts only small nuclear modification for mid-rapidity pions (see Fig. [3] below), while at 
RHIC energies it predicts large modifications, as CGC does, for forward produced pions. 



The much higher energy range at the LHC also opens the kinematic window on sma^ 
physics that can be explored in p + A collisions at mid-rapidity. Some CGC models [igI, 
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predict a suppression with Rpp\,{rj = 0,pt ^ lOGeV/c) ~ 0.5 with strong dependence on 
the initial evolution conditions. Such small values of Rpph would imply that nearly all 
nuclear suppression observed in NMF -RpbPb in Pb+Pb collisions, previously attributed to 
jet quenching in the final state, could instead be due to nonlinear initial state parton flux 
suppression. 
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Due to a factor of two increase in the final parton densities at the LHC, jet quenching 
is expected to produce higher suppression than at RHIC energies. Actually, the observed 
Pb+Pb suppression of pions at the LHC energy was surprisingly weaker than expected from 
RHIC constrained analysis extrapolated to the LHC [18|. Thus, from the perturbative final 
state interaction point of view there appears to be no room for initial state suppression. 
Therefore, a measurement of -Rppb at mid-rapidity significantly less than unity would con- 
tradict not only perturbative QCD (pQCD) models of the initial state nPDF evolution but 
also theory of the final state perturbative opacity series of jet energy loss. Since strong 
coupling AdS/CFT holography [l^ predicts even stronger final state suppression effects, 
an observation in p + Pb of significant deviations from unity would then call into question 
the validity of holographic interpretations of RHIC and LHC A + A results, including the 
applicability of minimal viscous hydrodynamics to apparent perfect fluidity. 

At sufficiently high energies and virtualities, QCD factorization theorems guarantee that 
jet observables can be calculated in perturbation theory. The open question is at what 
scale does factorization break down for nuclear processes. CGC theory 14l4l7l. I2CH28I] has a 
saturation natural scale Qs{x,A) that in principle provides the answer when Qs » Aqcd- 
However, nuclear jet observables up to LHC energies are sensitive to details of large x > 0.01 
as well as small A = 1 "corona" nucleon distributions for which ~ 1 GeV. 



Monte Carlo models as HIJINGl.O 



29j, HIJING2.0 |30j and HIJING/BB2.0 |3l| have been 



developed to study hadron productions in p + p, p + A and A + A collisions. They are es- 
sentially two-component models, which describe the production of hard parton jets and the 
soft interaction between nucleon remnants. The hard jets production is calculated employ- 
ing coUinear factorized multiple minijet within pQCD. A cut-off scale Pq in the transverse 
momentum of the final jet production has to be introduced below which [pt < Po) the 
interaction is considered nonperturbative and is characterized by a finite soft parton cross 
section cTgoft- Jet cross section, depend on the parton distribution functions (PDFs) that are 
parametrized from a global fit to data 30 1 . 

_Nucleons remnants interact via soft gluon exchanges described by the string models 
34| and constrained from lower energy e + e,e + p,p + p data. The produced hard 



33 



jet pairs and the two excited remnants are treated as independent strings, which fragments 
to resonances that decay to final hadrons. Longitudinal beam jet string fragmentations 
strongly depend on the values used for string tensions that control quark-anti-quark (qq) 
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and diquark-anti-diquark (qqqq) pair creation rates and strangeness suppression factors (7^). 
In the HIJINGl.O and HIJING2.0 models a constant (vacuum value) for the effective value 
of string tension is used, kq = 1.0 GeV/fm. At high initial energy density the novel nuclear 
physics is due to the possibility of multiple longitudinal flux tube overlapping leading to 
strong longitudinal color field (SCF) effects. Strong Color Field (SCF) effects are modeled 
in HIJING/BB2.0 by varying the effective string tensions value. SCF also modify the frag- 
mentation processes resulting in an increase of (strange)baryons which play an important 
role in the description of the baryon/meson anomaly. In order to describe p + p and central 
Pb + Pb collisions data at the LHC we have shown that an energy and mass dependence of 
the mean value of the string tension should be taken into account 3l| . Moreover, to better 
describe the baryon/meson anomaly seen in data a specific implementation of J J loops, has 



to be introduced. For a detailed discussion see Ref. [3l|. Similar result can be obtained 
by including extra diquark-antidiquark pair production channels from strong coherent fields 
formed in heavy ion collisions 32 1. 

All HIJING type models implement nuclear effects such as nuclear modification of the 
partons distribution functions, i.e., shadowing and jet quenching via a medium induced 
parton splitting process (collisional energy loss is neglected) 29|. In the HIJINGl.O and 
HIJING/BB2.0 models Duke-Owen (DO) parametrization of PDFs 351] is used to calculate 
the jet production cross section with px > pq. In both models using a constant cut-off 
Po = 2 GeV/ c and a soft parton cross section asoft = 54 mb fit the experimental p + p 
data. However, ioi A + A collisions in HIJING/BB2.0 model we introduced an energy and 
mass dependence of the cut-off parameter, po{s, A) 3l|] at RHIC and at the LHC energies, in 
order not to violate the geometrical limit for the total number of minijets per unit transverse 



area. In HIJING2.0 



30|] model that is also a modified version of HIJINGl.O 29|] the Gluck- 



Reya-Vogt (GRV) parametrization of PDFs 36| is implemented. The gluon distributions in 
this different parametrization are much higher than the DO parametrization at small x. In 
addition, an energy-dependent cut-off po('5) and o"so/t(s) are also assumed in order to better 
describe the Pb -|- Pb collisions data at the LHC. 

One of the main uncertainty in calculating charged particle multiplicity density in Pb -|- 
Pb collisions is the nuclear modification of parton distribution functions, especially gluon 
distributions at small x. In HIJING type models one assume that the parton distributions 
per nucleon in a nucleus (with atomic number A and charge number Z), fa/Ai^jQ"^), are 
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factorizable into parton distributions in a nucleon ( fa/N) and the parton(a) shadowing factor 

fa/Aix, Q"^) = Sa/Aix, Q'^)fa/N{.X, Q^) (1) 

The impact parameter dependence is implemented through the parameter Sq, 

5 6^ 



s,{h) = s,-\l-^] (2) 



where Ra = 1.12A^/^ is the nuclear radius. 

In HIJING/BB2.0 the shadowing factor for gluon and quark are assumed to be equal 
{SgiA{x,Q'^) = Sq/A{x,Q'^)) and are similar with those used in HIJINGl.O They were 
selected in order to fit the centrality dependence of the central charged particle multiplicity 
density at the LHC In contrast, in HIJING2.0 a much stronger impact parameter dependence 
of the gluon {sg = 0.22 — 0.23) and quark (sg = 0.1) shadowing factor is used in order to 
fit the LHC data. Due to this stronger gluon shadowing the jet quenching effect has to be 



neglected 30|. Note, all HIJING type models assume a scale-independent form of shadowing 
parametrization (fixed Q^). This approximation could breakdown at very large scale due 
to dominance of gluon emission dictated by the DGLAP \^ evolution equation. At Q 
= 2.0 and 4.3 GeV/c, which are typical scales for mini-jet production at RHIC and LHC 
respectively, it was shown that the gluon shadowing varies approximately by 13% in EPS09 
parametrizations 

III. COMPARATIVE STUDY OF MODEL PREDICTIONS 

Figure [T] shows HIJING/BB2.0 predictions of the global observables dNch/drj and RpPhiv) 
= {dN^^^)/dri)/{N^oiidN^^/dri) characteristics of minimum bias p + Pb collisions at 4.4^4 
TeV. The predictions for p + p are also shown. Minijet cutoff and string tension parameters 
Po = 3.1 GeV/ c and k = 2.9 GeV/fm for p + Pb are determined from fits to p + p and A + A 



systematics from RHIC to the LHC (see Ref. 311 for details). Note, these calculations 
assume no jet quenching. 

The absolute normalization of dNc^/ drj is however sensitive to the low pt < "2 GeV/ c non 



perturbative hadronization dynamics that is performed via LUND (33| string JETSET [34 1 
fragmentation as constrained from lower energy e + e, e+p,p+p data. The default HIJINGl.O 
parametrization of the fixed Ql = 2 GeV^ shadow function leads to substantial reduction 
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FIG. 1: (Color online) (a) HIJING/BB2.0 predictions of charged particles pseudorapidity distri- 
bution (dNch/drj) for minimum bias (MB) p+Pb collisions at AAA TeV. Solid YS curve includes 



fixed shadowing functions from HIJINGl.O 



29|, while the dashed NS curve has no shadowing. 



(b) Ratio RpPbiji) calculated assuming A'^coii(MB) = 6.4 

(solid histograms) of the global multiplicity at the LHC It is important to emphasize that 
the no shadowing results (dashed curves) are substantially reduced in HIJING/BB2.0 relative 



to no shadowing prediction with default HIJING/1.0 from Ref. [29|], because both the default 
minijet cut-off = 2 GeV/c and the default vacuum string tension = 1 GeV/fm (used 
in HIJINGl.O) are generalized to vary monotonically with centre of mass (cm) energy per 



nucleon y/s and atomic number, A. As discussed in [31[, systematics oi p + p and Pb+Pb 
multiparticle production from RHIC to the LHC are used to fix the energy {\fs) and the 
A dependence. Thus the cut-off parameter pq{s,A) = 0.416 ^/s^'^'^^ y^o.i28 QeV/c and the 
mean value of the string tension /t(s, A) = (s/sO)°'°^ y^o.ie? GeV/fm. The above formulae 
lead to po = 3.1 GeV/ c and k, = 2.9 GeV/fm at 4.4^4 TeV for p + Pb collisions. For p + p 
collisions at 4.4 TeV we use a constant cut-off parameter popp = 2 GeV/ c and a string tension 
value of Kpp = 2.7 GeV/fm. 

Note, even in the case of no shadowing shown in Fig. [H the increase to po = 3.1 GeV/ c 
from pq = 2 GeV/ c (value used in p + p at 4.4 TeV) causes a significant reduction by a factor 
of roughly two of the minijet cross section and hence final pion multiplicity. This reduction 
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of minijet production also is required to fit the low charged particle multiplicity growth in 
A + A collisions from RHIC to LHC (a factor of 2.2) [37 1. 
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FIG. 2: (Color online) (a) Minimum bias transverse momentum distributions at mid-pseudorapidity 
\rj\ < 0.8 predicted by HIJING/BB2.0 with (solid histogram) and without (dashed histogram) 



HIJINGl.O shadowing functions 



2^. The results for p + p collisions at 4.4 TeV (dotted histogram) 



are also included, (b) The mid-pseudorapidity nuclear modification factor of charged hadrons 
RpPh from HIJING/BB2.0 model. The solid and thin dashed histograms have the same mea ning 
as in part (a). They are compared to pQCD leading order (LO) predictions (dash dotted) 13] 



30l^nd to DGLAP evolved nPDF, EKS99 (dotted) I4|. 



IJ (KKT04) and CGC - rcBK model (thick sohd) from 



using HIJING2.0 shadowing functions 

Predictions of CGC model (thick dashed) 
□ 

Ref. [27|] are also included. 

We interpret this as additional phenomenological evidence for gluon saturation physics 
not encoded in leading twist shadow functions. The pt > 5 GeV/ c minijets tails are un- 
affected but the bulk low pr < 5 GeV/c multiplicity distribution is sensitive to this extra 
energy {^/s) and A dependence of the minijet shower suppression effect. It is difficult to 
relate po to saturation scale Qsat directly, because in HIJING hadronization proceeds through 
longitudinal field string fragmentation. The energy (i/i) and A dependence of the string 
tension value arises from strong color field (color rope) effects not considered in CGC plie- 
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nomenology that assumes factorized gluon fusion hadronization. HIJING 



ladronization of 



minijets is not via independent fragmentation functions as in PYTHIA [3J], but via string 
fragmentation with gluon minijets represented as kinks in the strings. The interplay be- 
tween longitudinal string fragmentation dynamics and minijets is a nonperturbative feature 
of HIJING type models. The approximate triangular (or trapezoidal) rapidity asymmetry 
seen in the ratio RpPhiv) sloping downwards from the nuclear beam fragmentation region at 
negative pseudorapidity t] < —5 toward l/Nco\\ in the proton fragmentation region [rj > 5) 
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39| and realized via string 



is a basic Glauber geometric effect first explained in Refs. 
fragmentation in HIJING. 

In Fig. [2] are displayed the predicted transverse spectra and nuclear modification factor 
for charged hadrons at mid-pseudorapidity, |?7| < 0.8. Including shadowing reduces -Rppb 
from unity to about 0.7 in the interesting 5 to 10 GeV/c region close to the prediction of 
Color Glass Condensate model UM (KKT04). A similar nuclear modification factor is found 
|13l | using leading order (LO) pQCD coUinear factorization with HIJING2.0 parameterization 
of shadowing functions 40|], GRV parton distribution functions (nPDF) from Ref. 36|], and 
hadron fragmentation functions from Ref. 41|. 

In stark contrast to the three curves near 0.7±0.1 from completely different dynami- 
cal modeling the standard DGLAP evolved global e+A fit nPDF (dotted curve labelled 
EKS99 ^) predicts near unity for transverse momenta above 5 GeV/c. The no shadowing 
HIJING/BB2.0 values (NS, thin dashed histogram) goes to unity above 5 GeV/c, but the 
nonperturbative string hadronization pulls the intercept at = near to 1/2 as constrained 
by the global triangular enhanced form of dN„ph/ dy relative to dNpp/dr] shown in Fig. [lb. 



Note, that the model BGK77 from Refs. 



38 



also predicts RpA{y,PT = 0) = 1 at the 
nuclear target rapidity and l/Nco\\{b) at the proton projectile rapidity. 

However, a recent new version of the CGC-rcBK model predicts essentially no shad- 
owing/saturation effects at ?7 = in contrast to both CGC-KKT04 [l^ and CGC-rcBK 



model from Ref. 16[. The absence of shadowing at midrapidity in the CGC-rcBK 
model is due to a phenomenological extra anomalous dimension 7, introduced to modify 
color dipole cross section crdipoieif') (r^)"^. This significantly steepens the pp transverse 
momentum distribution relative to the quadratic form (Tdipoie{f) used in CGC model 
(MV) 42] as required to reproduce LHC pp data. Recently, posible extra A dependence of 



this extra anomalous dimension has been proposed 28|]. It would be very surprising indeed 
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if future p + Pb data would show no evidence of shadowing with a Rppb ~ 1.0 at 77 = 
mid-pseudorapidity which could then be ascribed either to (i) rapid DGLAP evolution 
of shadowing in EKS09 6|] parametrization or to (ii) accidental cancellation of deep satura- 
tion effects due to an anomalous short distance behavior of the dipole cross section in CGC 
modeling. 
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FIG. 3: (Color online) Predictions updated at AAA TeV of Refs. Q results for central - 20% 
(b < 3.5 fm) p + Pb at mid-rapidity. The original predictions at 0.2^4 TeV for d + Au are also 



included. Compared are the results obtained with fixed shadowing functions HIJING2.0 



3C| 



with (b-dep) and without (b-indep) impact parameter dependence. Predictions with DGLAP 
evolved shadowing functions from Ref. 4| (EKS99) and Ref. [s] (EPS08) are also shown. The data 



are from PHENIX Collaboration 



43| 



What is significant about higher pt deviations from unity is that the nuclear modification 
factor in central A -|- A collisions is related to the minimum bias p + A hj simple Glauber 
geometric considerations, that can be expressed as 
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RAA{,y = 0,pT,b = 0) ^ {{RpAiy = 0,pr,min.bias)))^ (3) 

where the average is calculated over all impact parameters. Thus, -Rppb ~ 0.7 for mini- 
mum bias collisions implies a NMF -RpbPb ~ 0.5 in central Pb+Pb collisions before any final 
state interactions take place. The result for RpPhiPr) ~ 0.7 is similar with those reported 
recently with CGC type models [igI, (rcBK), albeit with huge error bars at mid-rapidity 
because of poorly known initial saturation conditions for p + p and near the surface of heavy 
nuclei. However, if the prediction of -Rppb ~ 0.5 turned out to be confirmed by the 
upcoming p + Pb measurements then, in central Pb+Pb collisions we expect a factor of 
roughly four suppression 0.25) in pions at transverse momenta of roughly 10 GeV/c. 
This fact would leave no room for final state interactions in matter 100 times denser than 
ground state nuclei. Needless to say this point alone underlines more the importance of 
measuring the px dependence of NMF {RpPhiPx)) at the LHC. 



In Figure|3]the updated predictions at AAA TeV [13| of RpPhiPr) in central(0— 20%;6 < 3.3 
fm) p + Pb collisions at mid-rapidity are shown. The message is similar to that obtained 
from Fig. [2l The standard coUinear evolved nPDF models from Ref. jj] (EKS99) and 
from Ref. p] (EPS08) predict only a slight deviation 10%) from unity, as discussed in 
detail in Refs. fliiol. Fixed Q"^ shadowing functions used in HIJINGl.O or HIJING2.0 models 
predict Rpphipx) = 0.6 ± 0.1 in the pt range 5 to 15 GeV/c, well below unity. Previous 



results at the RHIC energy 



13| for central (0 - 20%) d + An collisions at 0.2A TeV are 



also presented in comparison with PHENIX data j43|. At RHIC energy {0.2A TeV) all 
models predict approximately RdPhiPr) = 1 ±0.1. At this energy HIJING2.0 the shadowing 
is much weaker for pt > 5 GeV/c domain because this correspond to x > 0.05, which is 
over an order of magnitude larger than at the LHC energy. Taking an impact parameter 



30| for central p + Pb collisions results in 



dependence of shadowing function in HIJING2.0 
a further decrease of the NMF RpPbipr) by 15 - 20%. 

Finally, in Fig. H] we show the nuclear modification factor RpPb for inclusive charged 
hadrons {h^ + h~) at ^/s = 4.4 TeV obtained from different models for minimum bias (MB) 
(A'coii = 6.4) in Fig. |4|l and central (A'^cou = 12) p + Pb collisions in Fig. Hb at forward 
pseudorapidity (r/ = 6). For reference we include the NMF at mid-pseudorapidity |?7| < 0.8 
(solid histograms) predicted with HIJING/BBv2.0. Smaller suppression in the region of 
interest {5 < px < 10 GeV/c) than obtained within HIJING/BB v2.0 model is predicted using 
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FIG. 4: (Color online) (a) The minimum bias (MB) NMF of charged particles at forward pseu- 
dorapidity r] = 6, from HIJING/BB2.0 model (dashed histogram). The results are obtained with 



291 ]. They are compared to pQCD leading order 



shadowing functions from the HIJINGl.O model 

(LO) results at r] = 6(dash dotted) [l^ using impact parameter dependent (b-dep) HIJING2.0 



shadowing functions 



30( 1 and to predictions obtained with DGLAP evolved shadowing functions 



5|. For reference the results 



(dotted) with no impact parameter dependence (b-indep) from Ref. 
at mid-pseudorapidity, |r/| < 0.8 (solid histogram) are also included, (b) The results obtained for 
NMF of charged particles in central (A'coU = 12) p+Pb collisions at AAA TeV. The histograms and 
the lines have the same meaning as in part (a). 

EPS08 parametrization Q with no impact parameter dependence (b-indep) of shadowing 
functions (dotted) and with pQCD calculation using an impact parameter dependence (b- 
dep) from HIJING2.0 model. In this range of transverse momenta, the HIJING/BB v2.0 
model predict slightly higher values than predicted (0.35 -0.40) by CGC model (rcBK) 16 |. 

It is obvious that at forward pseudorapidity the suppression is higher for central than in 
minimum bias (MB) p + Pb collisions. Moreover, for central p + Pb collisions the sensitivity 
to parametrization of shadowing functions is amplified. The different shape predicted by HI- 
JING/BB v2.0 at forward pseudorapidity in both minimum bias (MB) and central collisions 
could be explained as a specific interplay between JJ loops and SCF effects embedded in the 
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model, which induce a baryon/meson anomaly. Note, the same effect has been predicted in 



p + p and Pb+Pb collisions at LHC energies [SjJ . To draw a definite conclusion, measure- 
ments of identified particle NMF, -RpPb(PT) are needed. Such measurements will provide a 
vital information on cold nuclear matter effects and will constrain the main parameters of 
shadowing functions used within different models. 



IV. CONCLUSION 

In conclusion, even with a small sample of 10^ events the study of RpPhiPx) or central 
relative to peripheral NMF (i?cp(pr)) could provide a definitive constraint on nuclear shad- 
owing implemented within different pQCD inspired models and CGC saturation models, 
with high impact on the interpretation or reinterpretation of the bulk and hard probes for 
nucleus-nucleus (Pb+Pb) collisions at LHC energies. 

HIJING type models predict for central (A'coii = 12) p + Pb collisions smaller values of 
nuclear modification factors (RpPhiPT)) than in minimum bias events. The possibility to 
trigger on the highest multiplicity tails of transverse momentum spectra in p + Pb collisions, 
will open the way to study collective phenomena in proton nucleus interactions with super- 
dense nuclear cores, where average number of binary collisions could increase to N^ow > 12. 
These measurements will provide an stringent test of the phenomenological models discussed 
in this paper. 
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